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^ It appears that at least several percent of large Kuiper belt objects are members of 
^ wide binaries. Physical collisions are too infrequent to account for their formation. Col- 
^lisionless gravitational interactions are more promising. These provide two channels for 
binary formation. In each, the initial step is the formation of a transient binary when 
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fS| two large bodies penetrate each other's Hill spheres. Stabilization of a transient binary 
^ requires that it lose energy. Either dynamical friction due to small bodies or the scat- 
^ tering of a third large body can be responsible. Our estimates favor the former, albeit 
^ by a small margin. We predict that most objects of size comparable to those currently 
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pbserved in the Kuiper belt are members of multiple systems. More specifically, we 
O . 

^ derive the probability that a large body is a member of a binary with semi-major axis 
^o,„.. X.e p..... „po„ ...ee „ ..ace 

^ density of large bodies having radius R, and 9q ^ 10 ^, the angle subtended by the solar 
?^ radius as seen from the Kuiper belt. For {o/Y?)R < a < R/9q, the probability is just 
{Tj/ pR)6q^, the optical depth of the large bodies divided by the solid angle subtended 
d by the Sun. For R < a < = {a/T.)R, it varies inversely with semi-major axis and 
reaches ~ {o / pR)6q^ at a ^ R. Based on current surveys of the Kuiper belt, we estimate 
S/p ~ 3 X 10""^ cm and R ~ 100 km. We obtain a/p 0.3 cm by extrapolating the surface 
density deduced for the minimum mass solar nebula. Rough predictions are: outside 
of the critical separation r^/a© ~ 3", the binary probability is ~ 0.3%; at separations of 
0.2", comparable to current resolving capabilities, it reaches ~ 5%, in agreement with 
results from the HST binary survey by Brown. 



2 mainauthor should go here 

1 Introduction 

The Kuiper belt^ is the best solar system laboratory for studies of the early stages of runaway 
accretion. Runaway accretion in the Kuiper belt terminated when the velocity dispersion of its 
members was increased by an as yet undetermined process.^ Unhke the asteroid belt, its largest 
members have suffered little coUisional evolution since. 

The discovery that a substantial fraction of its largest members are in binaries with wide sepa- 
rations and order unity mass ratios^"^ is the latest of many surprises provided by the Kuiper belt. 
Collisions coupled with tidal evolution, mechanisms that may explain other solar system binaries, 
fail to account for the formation of the Kuiper belt binaries. 

The low frequency of collisions among the large Kuiper belt objects^° implies that binaries formed 
by coUisionless interactions mediated by gravity. These were most effective earlier when dynamical 
friction due to small bodies limited the velocity dispersion of the large ones. This situation pertains 
during runaway accretion. 

In the following section §2, we outline a simple model for runaway accretion to set the stage for 
binary formation. We are guided by more detailed formalisms implemented in numerical simulations. 
Values of relevant parameters are estimated based on the numbers and sizes of objects deduced 
from Kuiper belt surveys, and by the extrapolation of the surface density in the minimum mass 
solar nebula. We estimate the binary formation rate and derive the semi-major axis distribution 
in section §3. In the final section §4, we discuss the observational implications of our results and 
mention several open issues. 

2 Preliminaries 

In this section, wc introduce a simple set of equations to describe the major processes that occur 
during the growth of planetesimals. These equations are abstracted from more complete treatments 
in the literature. ^^'^^ The main simplification in our approach is the identification of two groups of 
bodies, small ones, containing most of the total mass, and large ones, contributing a small fraction 
of it. The latter group, which dominates the stirring of all bodies, is identified with the currently 
observed population of Kuiper belt objects. 

Gravitational interactions determine the accretion rates and velocity dispersions of both large 
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and small bodies. The governing equations for these processes read: 
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Q^udt^ pR ^-"^ ps' ^ ' 

where R and s are the radii of the large and small bodies, E and u are their surface mass densities, 
V and u arc their velocity dispersions,* p is the material density^ and VLq is the orbital frequency 
around the Sun.-^- 

The expressions TD^q/ pR and (tVIq/ps give the kinematic (ignoring gravitational focusing) col- 
lision rates of large bodies onto large bodies and small bodies onto small bodies; gQ-q/ pR is the 
kinematic collision rate of small bodies onto a large one multiplied by the mass ratio (s/R)^. The 
factors Fm-m and FM-m in equation (1) account for the enhancements of the physical collision rates 
due to gravitational focusing by large bodies of the trajectories of incoming large and small bodies 
respectively. Cross sections for large deflections exceed those for physical coUisions when gravita- 
tional focusing is significant. The factors F^_j^ and Flj_^ that appear in equations (3) and (4) are 
the ratios of those cross sections to the geometrical ones. Appropriate expressions for Fm-tu read: 

{1 Vesc < U 

{Vesc/uf VH<U<Vesc (5) 
{'"Ischnu) Oq'^vh <u<vh 
where ^0 10~^ is the angle subtended by the solar radius as seen from the Kuiper belt, Vesc is the 
escape velocity from the large bodies and vh is their Hill velocity defined in terms of the Hill radius, 

*Wc are making an additional simplification by using a single velocity dispersion in place of three that are needed 
to characterize a triaxial velocity ellipsoid. 

'^We use the same density for the Kuiper belt objects as for the Sun. 
■tin writing equations (l)-(4) we implicitly assume that u> v. 
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as 

vh ~ ^qRh- (7) 
The escape velocity is related to the Hill velocity by 

Vesc-^iGpy/'R^VH/eH^ (8) 

Fm-m is obtained by replacing u with v in the above expressions and in their ranges of validity. 

Each of equations (l)-(4) has a simple interpretation. In equation (1), the two terms on the rhs 
account for the radius growth of the large bodies by the accretion of large and small bodies. The 
single term on the rhs of equation (2) describes how the radii of small bodies increase as the result 
of coalescence under conditions of negligible gravitational focusing. In equation (3), the first term 
describes the viscous stirring of large bodies due to mutual gravitational deflections, whereas the 
second term accounts for the damping of the large bodies' velocity dispersion by dynamical friction 
resulting from their gravitational interactions with the small bodies. Equation (4) shows how the 
velocity dispersion of the small bodies evolves under viscous stirring by large bodies and damping 
due to coUisions between small bodies. 

Values of the dimensionless parameters E/u <^ 1 and ^0 determine the appropriate regime of 
runaway accretion. At the Kuiper belt a© ^ 40 AU so 9q ^ 10~^. Based on Kuiper belt surveys^^"^^ 
which detect only the large bodies, we estimate S/p ~ 3 x 10^'^cm and R ~ 100 km. Extrapolating 
the minimum mass solar nebula surface density to 40 AU, yields a j p 0.3 cm. The resulting ratio 
of E/(7 ~ 10~^ is compatible with Kuiper belt simulations. ^^"^^ In the rest of this paper we assume 
that 

< - < 1. (9) 

(7 

In this regime of runaway accretion, the rates of gravitational stirring and dynamical friction 

acting on the large bodies are much greater than the large bodies' growth rates. Moreover, provided 
the radii of the small bodies satisfy 

s>-i?~ 0.1km, (10) 
(J 

their growth rates are negligible compared to those of the large bodies and coUisional damping of 

their velocities is unimportant on the timescale of viscous stirring.^ Under these approximations, 
§We assume that the inequaUty in equation (10) always applies. 
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dynamical friction on the large bodies balances their viscous stirring, and the stirring of the small 

bodies proceeds with negligible coUisional losses. Equations (l)-(4) yield 

1 /S\3/2 



V 



^0 

and 



~ 0.3 (11) 



-~(r-r-3- (12) 



Thus, the velocity ratio 



- - 0-,^'% ^ 10. (13) 

These results are in accord with the most recent simulations.^^ 

Given this solution, the balanced rates of viscous stirring and dynamical friction for large bodies 
satisfy 

"''^ '''' -^-^'Oe. (14) 



V dt 
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1 dv 
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Large bodies grow predominantly by accreting small bodies; The small bodies' velocity dispersion 
evolves mainly by viscous stirring provided by the large bodies. These two processes proceed at 
equal rates: 



1 dR 1 du (7/(7 
R dt u dt pR \E 



Oo- (15) 



3 Binaries 



Because v < Q,qRh, the gravity of a large body significantly deflects other large bodies at separations 
smaller than Rh, whereas small bodies are only affected at separations^ smaller than 

'r. (16) 



E 



Since vh < u < Vgsc, we have R < < Rh- 

For separations r < Vu, the velocity dispersion of the small bodies is 

C^f^Tf- (-) 

In addition, since the impact parameter for a small body to arrive at radius r is 6 ~ (r^r f\ the 
continuity equation implies that the small bodies' number density is enhanced by a factor 

h^n, /r \ ^/-^ 



s^w(r„/s)^/^ V s 

^We assume that the small bodies are on hyperbolic orbits with respect to the large body. 



0.03. (20) 
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3.1 binary formation 

We identify two distinct channels for binary formation. Both begin when two large bodies penetrate 
each other's Hill spheres. This occurs at a rate 

per large body. Stabilization of such a transient binary requires energy loss on timescale Oq^. This 
can be achieved either by dynamical friction from small bodies or by interaction with a third large 
body. 

Dynamical friction from small bodies during time Qq^ results in a fractional energy loss 

Mvjj ^ pR Ve; 

This is also the fraction of transient binaries that become bound. Numerical simulations that verify 
this assertion are described in the figure. Thus, this mechanism leads to a binary formation rate per 
large body 

The probability that a third large body joins the Hill sphere of a transient binary during its 
lifetime is 

A significant fraction of these triplets will result in a bound binary. Therefore, the binary formation 
rate per large body via this channel is 



With our parameters the ratio 



7^l Ka 

is smaller than unity. Thus from here on we take TZ = TZi. 
The ratio of 7^ to R'^dR/dt is 



(24) 



SO binaries form at a rate comparable to the large bodies' growth rate. 
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3.2 binary fraction and semi-major axis distribution 

Let p(a) be the differential probability distribution of finding a large body in a binary with semi- 
major axis a. In steady state" 

a p(a) ~ — , , , , , , (26) 

For < a < Rh^ a~^da/dt equals the energy decay rate due to dynamical friction given by equation 
(14). Thus 

which is independent of a. 

For R < a < two competing effects come into play. The number density of small bodies is 
enhanced by a factor (r„/ a)^/^ and their velocity dispersion is increased by the same factor. Therefore 
a~^|(ia/(it| is reduced by a factor a/r^ in comparison to its constant value for r„ < a < Rh- Thus in 
this interval of semi-major axis 

a ( iV R /R 



which increases inwards as a~^. 

The timescale for a binary to spiral in until contact is achieved is equal to R{dR/dt)~^. Contact 
occurs on the same time scale as that during which large bodies grow by accretion of small ones. 
Thus, if we define a critical separation 



dcrit ~ ~ 300 km, (29) 



inside of which a p{a) > 1, then accretion by binary inspiral might make a substantial contribution 
before R reaches Ucrit- As long as i? < Ucrit, we may also expect to find more than one companion 
to a given large body. However, both these assertions are uncertain since they require the stability 
of clusters of high multiplicity. For our parameters, R ~ acrit- 

Exchange reactions, in which the lighter member of a binary is replaced by a heavier body which 
passes through the system, are rare occurrences. The rate at which large bodies pass through an 
existing binary of semimajor axis th is given by equation (19). This is smaller than the rate of 
orbital decay. For smaller semimajor axes, even fewer large bodies pass through during an orbital 
decay time. 

II By steady state we mean during a time over which R is sensibly constant. 
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4 Discussion 



We propose that the wide binaries observed in the Kuiper belt formed during runaway accretion. 
A fraction of the large bodies that entered each other's Hill spheres became bound as the result 
of energy lost to small bodies by dynamical friction.** The time scale for a large body to become 
bound to a similar companion was 



Further dynamical friction hardened the binaries. The timescale to achieve contact was that 
during which isolated large bodies grew; 



Initially the inspiral was at a constant rate, but it slowed down inside r„ where the small bodies' 
number density and velocity dispersion were enhanced above their background values. We deduce 
that the probability that a large body is part of a binary with angular separation greater than 
ru/ttQ ~ 3" is 



Inward of r„, the binary probability per logarithmic interval of semimajor axis increases inversely 

with semimajor axis. Close to contact, the probability exceeds unity for R < Gcrit ~ 100 km, which 

implies that systems with higher multiplicities would exist if such are stable. For the resolution of the 

HST survey ~ 0.2", we predict a binary fraction of about 5%, roughly compatible with observations 

(Brown, private communication). Our prediction that close binaries are common could be tested by 

monitoring the brightness of Kuiper belt objects for evidence of echpses and/or fast rotation. The 

latter is a consequence, but not a unique signature, of binary mergers. 

In our analysis we estimate the surface density ratio of large to small bodies from the observational 

census of Kuiper belt objects together with an extrapolation of the surface density of the minimum 

mass solar nebula. It would be an improvement to have a theoretical understanding of how this 

ratio evolves during runaway accretion. Also, we elaborated the binary evolution scenario with a 

given set of parameters, specifically 9q and E/cr. Different scalings that apply for other parameter 

regimes remain to be worked out. 

**Interaction with a third large body was a less important channel for binary formation. 




(30) 




(31) 




(32) 
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Additional predictions suitable for testing by future observations of Kuiper belt binaries could 

be made by extending our model in a variety of ways. Relaxing the restriction to two mass groups 
would allow us to investigate the mass ratio distribution of binary components and its dependence on 
semi-major axis. We predict the formation of multiple systems for bodies with R < Ucrit- Numerical 
experiments can establish what types of systems survive disruption due to mutual interactions of 
their members. We could also predict eccentricity and inclination distributions. This would require a 
more detailed look at the anisotropy of the velocity distribution of small bodies, especially for a < r„. 
Anisotropy leads to tensorial dynamical friction in which the force is not antiparallel to the velocity. 
Binney^^ has investigated a related problem, the evolution of the orbits of galaxies in anisotropic 
clusters. We can already say that the eccentricities are likely to be large. They begin large following 
binary capture. Moreover, since the magnitude of the dynamical friction force increases outward for 
a < r„, they should remain large in this region. 

Simple numerical simulations of the Hill problem that incorporate dynamical friction would help 
pin down the dependence of the capture probability on the velocity dispersion and on the strength 
of dynamical friction. Computations involving three massive bodies arc needed for evaluation of 
capture probabilities via that channel. Three body capture cross sections have been computed for 
stars in globular clusters, but in our regime the Hill sphere must be taken into account. 

The two largest known Kuiper belt objects, Pluto and Charon, form a binary. CoUision coupled 
with tidal evolution is the generally accepted mode of formation for this system. Our scenario for 
binary formation provides another candidate. It may also apply to the formation of some binaries 
in the asteroid belt. 

In closing, we caution that our numerical results are crude. Dynamical processes have been 
estimated to order of magnitude. We have ruthlessly discarded numerical coefficients and logarithmic 
factors. More credence should be accorded to the functional dependences of our results than to the 
precise numerical predictions. 

Received 2 February 2008; Accepted draft. 
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Figure 1. Numerical simulations of binary formation by dynamical friction. Two equal mass bodies 
approach each other on initially circular orbits. We integrate the equations of motion including 
dynamical friction under the Hill approximation in a frame rotating at the average mean motion. 
The center of mass is fixed at (0,0). We plot only the trajectories of one of the bodies; the trajectories 
of the other body are related by reflection through the origin. Dashed red trajectories did not form 
binaries, and solid blue ones did. The axes are in 

unitsof a0(M/M0)i/3_ jj^g 

Lagrange points of the 

corresponding restricted 3-body problem are marked by green dots and the Hill sphere by a dashed 
green circle. Dynamical friction is modelled as a force antiparallel and proportional to the velocity 
measured relative to the Keplerian shear. Every capture orbit is plotted but only a sample of the 
non-capture orbits are shown. It is evident that capture is only possible within distinct ranges of 
impact parameters, whose widths are roughly proportional to the drag. The symbol C denotes the 
ratio of the linear measure of the impact parameters that result in capture to aQ^M/MQY^^. The 
symbol D denotes the fractional decrease in velocity due to drag over a time ^q^, essentially half 
the rhs of equation (20). Clearly C ~ D to better than a factor of two. This verifies our assertion 
that the fractional energy loss over time fl^^ is also the fraction of transient binaries that become 
bound. A sample of similar calculations starting from orbits with finite velocity v < vh gives similar 
results. 



